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ABSTRACT 



Current meter data were collected from moorings on the 800 and 1800 m isobaths 
on the continental slope off Point Sur, California at 100, 350, and 500 m depth and at 
1000 m depth on the 1800 m isobath. Complex demodulation and spectral techniques 
were used to study the distribution of tidal energy at the diurnal and semidiurnal fre- 
quencies, and to study the size, rotation and orientation of the tidal current ellipses. At 
mooring P2 (800 m) the observed diurnal energy increased with depth for each of the 17 
months of data analyzed, and the semidiurnal energy likewise increased with depth for 
14 of the 17 months. The observations were studied using the theory of bottom trapped 
waves in a rotating stratified fluid (Rhines, 1970). Bottom trapping is expected at the 
diurnal and semidiurnal frequencies whenever the critical trapping frequency 
= iVsinasin0 exceeds 1/24 and 1/12 hours'* respectively, where N = the Brunt-Vaisala 
frequency, a the bottom slope, and 0 is a measure of the angle the wavenumber vector 
makes with the gradient of topography. Theoretical energy decay curves matched the 
observed energy levels at P2 very well, indicating that the enhancement of energy with 
depth was likely due to the presence of bottom trapped waves. No bottom trapping was 
observed at P3 (1800 m) because the deepest current meters were positioned too far 
above the sea floor to detect bottom trapped energy, and because the density gradients 
near the sea floor were weaker than at site P2 and kept cu, below the tidal frequencies. 
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1. INTRODUCTION 



A. BACKGROUND 

The Naval Postgraduate School maintained an array of subsurface current meter 
moorings off Point Sur, California from 19S9 to 1990. The moorings were positioned 
on the continental slope and collected data used in the study of the California Current 
System (Figure I). The initial spectral analysis of data from mooring P2 indicated that 
the kinetic energy at diurnal tidal frequencies increased with depth (Figure 2). This was 
an unusual phenomena because vertical energy structures at tidal frequencies are usually 
associated with the baroclinic modes of free internal waves. At the latitude of Point Sur 
(36® 20' N) the inertial frequency is approximately 1.1 cpd and free internal waves can- 
not exist at diurnal (subinertial) frequencies. Baroclinic waves can exist at subinertial 
frequencies only if the waves are trapped along a wall or slope. 

The tides off most of the California coast are a mixture of diurnal and semidiurnal 
components. The principal lunar constituent, is the strongest tidal element followed 
by the and Sj constituents. Barotropic diurnal and semidiurnal tides both propa- 

gate northward along the coast. The semidiurnal tide is a combination of Kelvin, 
Poincare and forced waves. The barotropic semidiurnal tidal current ellipses, represent- 
ing horizontal current amplitudes and directions, are generally oriented parallel to the 
coastline {Munk ei al., 1970). The barotropic diurnal tide is a combined Kelvin and shelf 
wave with current ellipses which can be either parallel or perpendicular to topography 
{Noble el al., 1987). Several studies of California tides have found the primary direction 
of rotation for diurnal and semidiurnal current ellipses to be counterclockwise. 

Baroclinic tides are generated when semidiurnal surface tides force isopycnal sur- 
faces over topographic features, such as the continental shelf break {Wunsh, 1975). 
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Wliilc the barotropic semidiurnal tide is usually aligned with topography, internal tides 
often become the dominant cross-slope current component along the California coast 
{Noble cl ai, 1987). 

Variations in the vertical distribution of tidal energy have been observed at 
semidiurnal frequencies. During the Coastal Ocean Dynamics E.xperiment (CODE), the 
baroclinic semidiurnal internal tide was larger than the barotropic tide. Changes in 
semidiurnal kinetic energy at different depths were linked to surface winds and changing 
stratification {Rosaifcld, 1990). The semidiurnal internal tides observed in CODE had 
some bottom energy intensification which was attributed to the horizontal velocity 
structure of the first dynamic mode. 
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There are few published studies which discuss similar bottom intensification of 
subinertial currents. Bottom trapped topographic Rossby waves were observed by 
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Thompson and Luyten (1976) and were discussed using bottom trapped wave theory 
developed by Rhines (1970). These were low frequency waves, however similar theory 
was applied to bottom trapped waves at tidal frequencies and used to study the tidal 
currents and vertical energy distribution off Point Sur. 

B. BOTTOM TRAPPED WAVES IN A ROTATING STRATIFIED FLUID 

This discussion follows the development presented by Rhines (1970) for bottom 
trapped waves in a uniformly stratified, rotating fluid bounded by a single plane wall 
inclined everywhere at angle a to the horizontal. A right handed orthonormal coordi- 
nate system is oriented on the slope such that the x^ and jtj axis are rotated at an angle 
4> from the alongslope and upslope directions respectively and is perpendicular to the 
bottom slope (Figure 3). The corresponding velocities along Xi_ atj and Xj are designated 
as w, , and Uj . 

The governing equations are the linearized, non-diffusive Boussinesq equations used 
in the form: 



du , 

-57 + 22u = --^-«PJ 



(I) 



op 

di 



g 



uj 



( 2 ) 



where p = pressure, Q. = angular rotation of the earth, g = gravitational acceleration, 

A A 

j = a vertical unit vector, 2£2 = yj = /(sin a sin 0, sin a cos 0, cos a) , /= 

—g dpi 

Coriolis parameter, the Brunt-Vaisala frequency squared = PP = — ^ — , and a per- 

Po 01 

turbation density is used such that p = Po + P{^) + p'{x, y, z, t) . 

Inviscid boundary conditions are chosen so that Uj = 0 at Xj*0. every- 

where in the water column the motion is nondivergent in planes perpendicular to the x, 

5 «, du2 . . . . 

axis and the continuity equation -z 1- — — = 0 applies. The wave propagates m the 

0 X 2 
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Figure 3. Coordinates for Bottom Trapped Waves on a Sloping Bottom: The ele- 

vation from the horizontal is given by a. 



^u^ du^ 

Xj direction with — = 0 . To balance the continuity equation — — must also equal 

0X2 

d d 

zero. When u, # 0 and — = 0, free internal waves can exist. If Uj = 0 and -r — ^ 0 

0x2 0x2 

waves are trapped against the w'all (bottom) in a plane parallel to the boundary surface 
wdth particle motion limited to the x, direction. 

Substituting Uj, t/j = t /3 = 0, and the appropriate form of 2Q into Equations (1) and 
(2) leads to a set of three equations where for the x, direction, Equation (1) becomes 
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( 3 ) 



OM, 

— — = —gp'( sin a sin 0) 
cl 



In the X 2 direction. Equation (1) becomes 



, N \ f . Sp . op 

/«,( cos «) = - — ( ^ + -^ + 



Po \ o.r, 8 x 2 dxj 



-gPl 



1 op 

/w,( cos a) = - — - gp'{ sin a cos 4 >) 



(4) 



For the Xj direction, the result is 



1 op 

-/«,( sin a cos 0) = - — - gp'{ cos a) 
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(5) 



The rate of change of the density perturbation, equation (2), similarly reduces to 



Sp' _ 
dr S 



t/j( sin a sin 0). 



( 6 ) 



The combination of equations (3) and (6) results in an expression for the frequency a 
where: 



dp' 

8 i 



du^ 

6t 



= {N^Uip'){ sin a sin 0)^ 



{p'uy) = sin a sin 4 >)^ 

8t 



o' = sin a sin </>)^ 
<7 = + A* sin a sin 4 > 



(7) 
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Substituting (7) into (6), p' can be deternained as 



, .'V 

P = — u, 

Solutions to the wave equations are found using the form 



(8) 



( 9 ) 



where k and are real, k is defined as a trapping coelTicient with units of wavenumber. 
The alongslope or horizontal ( ) and vertical ( k -^ ) wavenumbers in this solution lead 

to spatial dependence in the and directions. Using equations (8) and (9) in (4) and 
(5), the spatial dependence of the wave is: 




SCl — sin^a sin^0] ^ 

5 -:r ~2 TT sin a sin 4>) 

cos a + 5 sin a cos (p 



( 10 ) 



k^ 1 (1 — S^) sin 2a cos 4> 

2 cos^a + sin^a cos^4> 



N 

where S = — . When the trapping coeflicient k is positive the wave motion is expo- 
nential and trapped near the bottom. From equation (9), wave motion decreases due to 
the effect of trapping by 



( 12 ) 

The magnitude of trapping is dependent upon a and <f). When particle motion is di- 
rectly up the slope (0 = 90°) trapping is maximized. When particle motion is along the 
isobaths (0 = 0°), no trapping occurs. 

Rhines discusses four difierent geometries where trapping occurs: (1) a ~ 1 for steep 
slopes, (2) a ~ S~' for more gradual slopes such as those on the continental slope, (3) 
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a S-' for very gradual slopes such as those on the continental rise, and (4) 
— a < S cos </) for near vertical slopes. Consider the ease appropriate for the conti- 
nental slope olT the California coast where a ~ 5'' and sin 0(1). Evaluating 
equation (7) using these relationships leads to the following; 



a 



± .V sin a sin (j) ~ A’ sin(S ') sin 0 




o f 



( 13 ) 



— can be expressed as: 

K-y 



K- -SCl — sin^5 * sin^(/)3 

^2 cos^S~' -I- sinks'"' cos^</> 



= 5 



( 14 ) 



The ratio of reduces to 



Aj 5Cl — sin^(5 ) sin^(/>D 

^2 cos‘(S~') -f- siii^(S“*) cos^0 



= S 



(15) 



The relationships now established with a ~ 5"' arc a ~ ~ S, and 

A, . . '. . 

— — - S. Parallel to the slope Coriolis is balanced by buoyancy as shown in equation 

(15). Coriolis traps waves against the wall (bottom) but the frequency of these trapped 
waves, determined by equation (7), is not dependent on /. Bottom trapped internal 
waves arc similar to internal Kelvin waves because they are not bound by a lower fre- 
quency limit of /and trapped waves can exist and propagate at subinertial frequencies. 
The high frequency limit of bottom trapped waves is a and is a function of a and i\. 
Maximum trapping will always occur when <t> = 90° and the particle motion is directly 
up the slope. 



8 



The exponential decay of trapped waves away from the bottom is determined by 
trapping cocfllcients. To determine k, bottom trapped waves were assumed to have a 
dispersion relationship similar to baroclinic Kelvin waves (Rhines, 1970) where 



with CO the frequency of a wave, k a horizontal wave number, n the mode number and 
H the water depth. The use of this dispersion relationship represents an approximation 
because baroclinic Kelvin waves assume a flat ocean bottom. The horizontal wave 
number can then determined using 



When 4> = 90° waves are propagating along the slope in the Xj direction and particle 
motion is up and down the slope in the jt, direction. Trapping coeflicients are defined 
as K = kjS from equation (14) and can be calculated using 



Waves are trapped closer to the sea floor and decay more rapidly aw'ay from the 
bottom when k increases. Trapping coefficients increase when (1) the buoyancy fre- 
quency N increases, (2) depth H decreases, (3) higher frequencies co are considered (in 
the tidal range of 1.0 to 2.0 cpd), (4) f decreases (lower latitudes) and (5) for higher 
modes n. Trapping is weakest and waves decay more gradually aw^ay from the bottom 
in the first baroclinic mode. No trapping occurs for the barotropic case when n = 0. 




rt= 1, 2, 3 ... 





(16) 
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II. DATA COLLECTION 



A. CURRENT METER DATA 

Current data were collected at two sites oil' Point Sur, California (Table 1) located 
over the continental slope (a = .06). 



Table 1. POINT SUR MOORING LOCATIONS 



Mooring 


Position 


Depth 


P2 


36-20.0 N 122-10.2 W 


800 m 


P3 


36-20.0 N 122-27.6 W 


1800 m 



Mooring P2, on the 800 m isobath, was equipped with three Aanderaa vector averaging 
current meters (RCM8) at depths of 100 m, 350 m and 500 m. Mooring P3 was situated 
52 km ofl'shore in 1800 m with four Aanderaa current meters at depths of 100 m, 350 
m. 500 m and 1000 ni (Figure 4). 

Vector averaged current velocities, current directions, temperature, conductivity and 
pressure measurements were recorded ever>' 30 minutes in a solid state Data Storage 
Unit (DSU). The DSU had a storage capacity of 65530 10-bit words and could record 
seven months of data with the sampling interval used. Sampling times were regulated 
by an internal quartz clock with an accuracy of ±2 seconds per day. Several deploy- 
ments were necessitated by the seven month data storage limit (Figure 5). P2 was re- 
covered and redeployed in May, August and December 1989 and May 1990. The three 
current meters used on P2 provided nearly continuous data from May 1989 through 
October 1990. Mooring P3 was deployed in December 1989 and again in May 1990. 
Good data was obtained at four levels from December 1989 to April 1990. After April, 
the data from the current meters at 500 m and 1000 m was limited by instrument failures. 



10 




Figure 4. Current Meter Configuration at Moorings P2 and P3 



Velocities were measured using a shrouded paddlewheel at the top of the recording 
unit. The paddlewheel was magnetically coupled to an electronic counter and velocities, 
in units of cm see"', were determined by the number of revolutions the paddlewheel 
made during each sampling interval. The sensor had a range of 2 - 250 cm sec*' with 
an accuracy of ± 1.0 cm sec"' or ± 2 percent of the current velocity. The paddlewheel 
stalled at velocities below a threshold of 2.0 cm sec"' (Aancleraa, 1990). The cosine re- 
sponse was acceptable within ± 27° from the vertical. Both the pressure sensors and the 
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POINT SUR CURRENT METER DEPLOYMENTS 


MOORING 




P2 100M 




350M 




500M 




P3 100M 




350M 




500M 




1000M 




JRN FZB ilflR 


RPR .'IRY JUN JUL RUG SlP CCT NOV OEC ^‘RN FE3 NRR RPR nRY JUN JUL RUG SEP OCT 


1909 


1990 



Figure 5. Current Meter Deployments at Moorings P2 and P3: The shaded bars 
indicate the intervals where good data was obtained at each depth. 



mooring dynamics program in use at NFS indicated a maximum mooring tilt of 11.7® 
which was well within acceptable limits. 

The current meters were oriented into the current by a large vane with an area of 
2425 cm^ The direction sensor was a magnetic compass with a potentiometer located 
at the bottom of the recording unit. The compass was oil damped and required 3 to 5 
seconds to resolve changes of 90° or more. The resolution of the compass was 0.35° 
with an accuracy of ± 5° for velocities 5 - 100 cm sec"' and + 7.5 degrees for velocities 
of 2.5 - 5.0 cm sec"' and 100 - 200 cm sec"'. The compass could function properly with 
a maximum tilt angle of 12° from the vertical. Prior to each deployment the combined 
magnetic deviation and variation of each unit was determined on a surveyed test bench 
at the Naval Postgraduate School and the necessary corrections were applied via 
postprocessing software. 
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B. HYDROGRAPHIC DATA 



While the current meters were deployed, several oceanographic research cruises were 
conducted in the region by the Naval Postgraduate School. These voyages consisted of 
California Undercurrent (CU) and Student Cruises (ST) which collected CTD data at 
intervals of approximately 10 km along the Point Sur Transect (POST). Data were col- 
lected near the moorings from the surface to within 50 - 150 m of the bottom using a 
Neil Brown Instruments System (NBIS) Mark IIIB CTD and are described in detail in 
Tisch (1990). Using the data from the cruises (Table 2), the various characteristics of 
the waters off Point Sur were examined. Although the station locations and spacing 
varied slightly with each cruise, the CTD data necessary for interpretation of the current 
meter results were comparable. 



Table 2. HYDROGRAPHIC CRUISES OCCURRING DURING 
CURRENT METER DEPLOYMENTS 



Cruise Designation 


Date of CTD Casts 


STMAY89 


5; 5 '89 


CLJUL89 


7;' 29, 89 


CUSEP89 


9 '26; 89 


STNOV89 


11'2 89 


CUJAN90 


1/18/90 


CUMAR90 


3/2 '90 


STMAY90 


5/2,90 


CLJUN90 


6'27'90 
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III. METHODS OF ANALYSIS 



A. SMOOTHING AND FILTERING 

Upon recover\’, the current meter DSUs were downloaded into a computer. The raw 
data files were converted standard scientific units using the sensor calibration equations 
and stored as formatted ASCII files. All data records were plotted to facilitate visual 
inspections. Velocity records were truncated where instrument failures or malfunctions 
were evident. Spurious values, such as an occasional outlier, were edited or deleted 
manually. 

All data were band pass filtered prior to further analysis to isolate the frequencies 
of interest. The filtering and interpolation methods used were originally developed by 
Oregon State University {Denbo ei ai, 1984). Postprocessed records, containing data 
at 30 minute intervals, were first low passed filtered (LP) using a Cosine-Lanczos filter. 
Such a filter minimizes unwanted ringing or oscillations caused by the Gibbs phenomena 
while still providing a sharp roll off at the desired cutoff frequency. A centered 25 point 
segment of the records was used to produce each LP filtered data point. The LP filter 
had a half power period of 2.9 hours (or a frequency of 8.4 cpd). The filter allowed any 
signal with a period exceeding 3.6 hours to pass through unattenuated. Signals with a 
period of less than 2.0 hours were removed entirely. After the LP filtering, records w'ere 
interpolated to transform 30 minute data records to 60 minute intervals with each data 
point falling on an even hour. The interpolation routine utilized Lagrangian 
polynomials to interpolate any unevenly spaced data points (due to missing individual 
data points or outliers) to specific time intervals {Gerald and IVhcaily, 1989). After in- 
terpolation, the data w'as again low passed filtered (LLP). The second Cosine-Lanczos 
filter had a half power period of 46.59 hours (or a frequency of 0.515 cpd) and used a 



14 



centered 121 point segment to produce one filtered data point. This filter was designed 
to remove diurnal and shorter period tidal energy and any inertial frequency energy 
present in the signal. The LLP records were subtracted from the LP records to produce 
the band passed form used in the tidal analysis. A sample from 100 m at P2 (Figure 6) 
shows the results of the filtering procedures. 

B. TIDAL COMPONENTS 
1. Harmonic Analysis 

To understand the nature of the tidal currents, an analysis of the various tidal 
components was desired. Traditional harmonic tidal analysis, or the technique of linear 
combinations, requires a complete set of observations over a standard length period. 
Tidal heights or motions can be expressed as a sum of harmonic terms, each determined 
by a known, fixed frequency. A single estimate of the amplitude and phase are then 
extracted through the analysis of observations. Observational periods are normally 
multiples of synodic periods involving the most important tidal constituents. The pre- 
ferred record length for a complete tidal analysis is 369 days. This allows extraction of 
the amplitude and phase of at least 64 of the largest tidal constituents {Dronkers, 1964). 
An adequate description of the major constituents can be accomplished using a much 
shorter observation of 29 days. A 14 day time series is adequate to provide the resol- 
ution of the principal diurnal constituents, AT, and 0^, and a 15 day time series will re- 
solve the principal semidiurnal constituents and 5^. Some of the longest tidal 
constituents require a full 19 year record for accurate resolution (Schureman, 1958). 

The current records collected off Point Sur var>' in length and are not ideally 
suited for the method of linear combinations. The traditional method could provide only 
a few independent estimates of the tidal constituents for each month of recorded current 
observations. Complex demodulation, a least squares analysis technique, was used in- 
stead to accomplish the harmonic analysis. The primary' advantage realized with this 
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CURRENT METER VELOCITIES P2 lOCM F12S9 




Figure 6. Current Meter Record Filtering: Time series from mooring P2 at 100 

m showing velocities in cm sec"* with LP filtering, LLP filtering and the 
band passed residual for the u (top three panels) and v (lower three 
panels) components. 
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method is that records of any length can be analyzed. Complex demodulation utilizes 
a set of normal equations which arc solved for unknown amplitudes and phase corre- 
sponding to n frequency constituents. The number of unknowns solved for in the set 
of equations is 2n + 1. The number of equations used in each set is determined by the 
number of observations, m, utilized for each estimate and the solutions are determined 
using matrix algebra. Schematically, this process an be represented as 

[observed velocities] = [known tidal constants][unknown empirical constants] 

The least squares method minimizes the difference between an observed current velocity 
and the sum of the solutions corresponding to that observation. The difference, or error, 
can be attributed to unresolved tidal constituents, inertial currents or noise in the signal. 
The extracted solutions reflect only those current variations with a coherent phase and 
no assumptions are made about data outside the interval upon which the solutions are 
fit (Pugh, 1986). 

The analysis of diurnal and semidiurnal constituents was conducted starting 
every six hours over 48 hour segments and produced a large number of independent es- 
timates for each current record. The 48 hour period did not allow adequate separation 
of the O, and A] diurnal and the A/j and Sj semidiurnal constituents which had prevented 
by synodic periods of 13.661 and 14.765 days respectively. Therefore the major tidal 
constituents (O,, A,, Q^, A/„ 7,; M 2 , S 2 , Aj, and A]) were generalized as diurnal ( 

o> ~ 1 cycle per day) or semidiurnal (o> ~ 2 cpd). Longer period constituents, such as the 
Lunar Long Period Mm and MSf, were removed by low pass filtering. Harmonics such 
as the Lunar Terdiurnal or Solar Terdiurnal S, were not considered due to the rela- 
tively small amplitude of these constituents w'hen compared to the primary diurnal and 
semidiurnal constituents. Higher order harmonics were removed by the initial band pass 
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filtering. Inclusion of the inertial frequency (o) = 1.1 cpd) as a known frequency pro- 
vided a better fit of the least squares solutions and reduced error. 

The complex demodulation program calculated the coefficients using a least 
squares approximation to the set of m equations with 



where (o is the angular frequency of the tidal constituent, i is the time of each hourly 
record and /, is the time of the first data point of an entire current meter record. The 
number of terms can be expanded to to include all desired constituents n. For u 
velocities, data was placed in a single row matrix U . The coefficients of equation (17) 
are placed in matrix A' . The following relationship was used to obtain the solutions 
to the set of equations represented by equation (17): 



Similar calculations were done for the v velocities to create matrix A V. Once the coef- 
ficients A were determined (the values in matrices AU or AV), the amplitude and phase 
of each frequency component was calculated using 



u OT V = A^ + A 2 sin (o{i — rj -b A^ cos cu(r — /^)... 



(17) 



(b'x A^‘) X (A'x A^') = /lf/ 





With these values, the current velocities can now be defined as 



M or V = Amp„ sin(cu(/ — + 6) 
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Diurnal and semidiurnal amplitude and phase were extracted from each current 
meter record using the complex demodulation technique. An example is provided in 
Figure 7. 

2. Construction of Current Ellipses. 

The amplitude and phase of clockwise (<2c, dc) anticlockwise 0„) ro- 
tating vectors were computed using monthly averages of the component amplitude and 
phase vectors (U,6J and (V,9,) obtained from complex demodulation using standard 
equations {Pugh, 1986): 



2c = y +V^~ 2UV sin(0, - 0j]T 
Qcc = + 2UV sin(0, - 



O^. = arctan< 



U sin 6^ + Fcos 1 
V cos — V sin J 



= arctan 



I 



V sin + F cos 9y | 
— U cos Og — V sin dy j 



Using the rotary components, the dimensions of a current ellipse may be defined as 

Semi-major Axis = <2c + Qac 
Semi-minor Axis = I <2c ~ Qac I 

Orientation of Semi-major Axis (clockwise from North) = -^ ^ {6^^ -I- 6^) 

Current rotation is clockwise if Q,. > Q„ and anticlockwise if > Qc- Sample ellipses 
are shown for mooring P2 during December 1989 (Figure 8). 
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AMPLITUDE AND PHASE: P2 500M PI 289 K1 




ure 7. Complex Demodulation of the Diurnal Frequency at 500 m on Mooring 
P2: The first two time series plotted are amplitudes in cm sec"‘ and the 
next two curves are the phases of the u and v diurnal velocity compo- 
nents in degrees relative to the beginning of each record. 
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TIDAL CURRENT ELLIPSE 

MOORING P2 DECEMBER 1989 




Figure 8. Tidal Current Ellipses: Current ellipses are shown for diurnal and 

semidiurnal tidal currents at P2 for December 1989. North is towards 
the top of the page. 
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Current ellipses from different moorings and varying depths were evaluated by 
visual comparison of the ellipse orientation, size, eccentricity and rotation and by 
quantitative comparison of the energy represented by each ellipse. Energy was com- 
puted in units of cm’ sec"’ by the addition of the squared dimensions of the semi-major 
and semi-minor axis. 

C. SPECTRAL ENERGY ANALYSIS 

To evaluate the energy present in the tidal currents, a program which computes the 
spectral characteristics of a time series was employed. The program, VCSPC3, was used 
to calculate the kinetic energy (KE) spectra as the total KE spectra for both the u and 
V velocities of a current meter record over one month periods. The program utilized a 
Cooley-Tukey fast fourier transform (EFT). One month records w'ere divided into five 
sequential, equal length segments. A separate spectrum was computed for each segment 
and the five segments were then averaged to produce the total energy spectrum. A 
parabolic data window was used in the program to eliminate discontinuities at the be- 
ginning and end of records created because record lengths were not exact multiples of 
the periods of each wave present in the signal {Bendai and Piersol. 1986). The effect of 
the data window was to improve the resolution of the narrow spectral peaks in the tidal 
frequencies by reducing spectral side lobe leakage and associated distortion. No addi- 
tional smoothing of the spectral estimates was done. The program calculated 95% 
confidence intervals using a Chi squared distribution. 

Program VCSPC3 provided both graphical output and data listings. The plots show 
kinetic energy density in units of cm’ sec"’ cpd*' versus frequency in cpd. Confidence 
intervals are indicated by the separation of parallel horizontal lines above the horizontal 
axis. Distinct peaks were clearly seen at both the diurnal and semidiurnal frequencies 
of these plots (Figure 9). 
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Figure 9. Spectral Energy at P2: A representative plot from the spectral analysis 
program is shown for a one month interval. The 95% confidence inter- 
vals is plotted below the curve. 
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The listings show spectral density at the diurnal and semidiurnal frequencies with a 
frequency resolution of 0.05 cpd and the associated 95 percent confidence intervals. 

D. HYDROGRAPHIC DATA 

A series of temperature (7), salinity (5), density (a,), geostrophic height (O), 
geostrophic velocity (V), and spiciness (;r), values were plotted in contoured vertical 
sections which passed through the locations of moorings P2 and P3 (figures not shown). 
A level of no motion for the geostrophic velocity calculations was chosen at 1000 m 
except where bottom depths were less than 1000 m when the greatest common depth 
betw^een stations was used. Brunt-Vaisala (or buoyancy) frequency N was computed for 
each cruise at the stations nearest moorings P2 and P3. Usually the CTD cast was 
within a few kilometers of the moorings' assigned position. Values for N were deter- 
mined by 



N 



2 




1 Sp 
P dz 



(18) 



Values of p used in equation (18) were computed using the International Equation of 
State 1980 for seawater {Millero and Poisson, 1981) using S, T and P values from the 
CTD cast and were smoothed by center averaging over 10 m increments. Calculations 
of used Az = 10 m. The results from each CTD cast were plotted as a vertical profile 
of A'^ against depth. Most A^^ profiles revealed the main pycnocline at depths between 
50 m and 200 m. The highest value of A’, 2.147 x lO'^cps (a 7.76 minute period), was 
found near the main pycnocline. Minimum N values w'ere found at the bottom of each 
cast with values near 2.344 x 10-^cps (a 71.01 minute period) and were assumed to be 
constant from the deepest calculated value to the sea floor. 

Profiles of the maximum trapping frequency tu,. = A’sina were constructed with the 
trapping frequency plotted in units of cpd against depth in meters. The possibility of 
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bottom trapping at any particular CTD station is best indicated by the magnitude of the 
maximum trapping frequency calculated over the lowest 10 m of CTD data (Figure 10). 
The deepest calculated typically 50 m off the bottom, was extrapolated to the bottom 
using the assumption that remained constant. Any wave with a frequency less than 
or equal to o, may be trapped. When is greater than 2.0 cpd, trapping may occur 
at both diurnal and semidiurnal frequencies. When is less than 1.0 cpd, no trapping 
of tidal currents should be expected. 
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Figure 10. Brunt-Vaisala Frequency and Trapping Frequency Profiles: 

and /V sin a are plotted versus depth using hydrographic data ob- 
tained near P3 in March 1990. 
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IV. RESULTS 



A. CURRENT ELLIPSES 

Monthly averaged tidal current ellipses were constructed using the amplitudes and 
phases obtained from the complex demodulation of current meter records (Figures 11 
to 14). The plotted ellipses display a great deal of variability in their size, orientation 
and direction of rotation. The size and energy represented by the current ellipses was 
found to be greatest for both tidal frequencies at the deepest current meter of mooring 
P2. 

I. Diurnal Current Ellipses 

Diurnal current ellipses for mooring P2 were constructed for seventeen consec- 
utive months (Figures II and 12). The energy represented by the ellipses was calculated 
using the sum of the squared axis lengths to compute kinetic energy in units of 
cm^ sec'^ . At the 100 m level the average energy of the ellipses was 18.0 cm^ sec~^. The 
direction of rotation of 76% of these ellipses was counterclockwise. At 350 m, the av- 
erage energy dropped to 14.5 cm^ sec”^ with 70% of the ellipses rotating counterclock- 
wise. At 500 m, the average ellipse energy increased to 22.3 cm^ sec"^ with 65% rotating 
counterclockwise. The diurnal current ellipses were larger and more energetic at 500 m 
than at 350 m for each of the seventeen months analyzed and there were twelve months 
when more energy was present at 500 m than at 100 m. During July 1989 the diurnal 
current rotation at all three P2 current meters was clockwise and there were eight 
months when the rotation at three depths was uniformly counterclockwise. 

The current meter records for six months were used to construct diurnal current 
ellipses at mooring P3 (Figure 12). At the 100 m level the average energy represented 
by the ellipses was 25.1 cm^ sec~^ with 66% of the ellipses rotating counterclockwise. 
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Figure 1 1. Diurnal Current Ellipses from Ma^' 1989 to November 1989 at Mooring 
P2: Horizontal ellipses are shown with north up and the direction of 

rotation indicated by arrows. A velocity scale is indicated at the bot- 
tom. 
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DIURNAL CURRENT ELLIPSES 



Figure 1 
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Diurnal Current Ellipses from December 1989 to September 1990 at 
Moorings P2 and P3: Details are as in Figure 11. 
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Figure 13. Semidiurnal Current Ellipses from May 1989 to November 1989 at 
Moorings P2 and P3: Details are as in Figure 1 1. 
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Figure 14. Semidiurnal Current Ellipses from December 1989 to September 1990 
at Mooring P2 and P3: Details are as in Figure 11. 
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V ORIENTATED NORTH 



At 350 m, the average energy dropped to 12.4 cm^ sec’^ . At 500 m, four months were 
analyzed and resulted in an average energy of 1 3.5 cm^ sec'^ with 75% rotating 
counterclockwise. The smallest ellipses and lowest energy (9.5 cm^ sec"^) were found at 
1000 m with 66% of the ellipses rotating counterclockwise. The ellipses were larger and 
contained more energy at 500 m than did those at either 350 m or 1000 m. There were 
four months when the current rotation at each P3 current meter was counterclockwise. 

The bathymetric contours at mooring P2 are oriented at approximately 005 ®T 
and the average orientation of the 500 m diurnal tidal current ellipses at P2 was 049° 
T. The semi-major axis of the diurnal ellipses crossed the isobaths a mean relative angle 
of 44°. At mooring P3 the bathymetry is oriented at 020° T. The average orientation 
of diurnal tidal ellipses at 1000 m was 074° T and the semi-major axis of these ellipses 
crossed the isobaths at a mean relative angle of 54°. 

2. Semidiurnal Current Ellipses. 

The semidiurnal current ellipses at P2 and P3 were larger and more energetic 
than those at diurnal frequencies. The orientation of the semidiurnal ellipses was less 
variable than the diurnal ellipses with almost all of them oriented in a north-northeast 
direction. 

At mooring P2, the seventeen month average energy represented by the ellipses 
at 100 m was 59.1 cm^ sec'^ The rotation of 58% of these ellipses was counterclockwise. 
At 350 m, the average energy decreased to 48.2 cm^ sec~^ and 52% of the currents ro- 
tated counterclockwise. At 500 m the average energy increased to 85.4 cm^ sec"^ with 
47% of the ellipses displaying counterclockwise rotation. The semidiurnal currents dis- 
played the same trends exhibited by the diurnal current ellipses where the largest ellipses 
and greatest energies were found at the deepest level and the smallest ellipses and least 
energy was at 350 m. The rotation of semidiurnal current ellipses at P2 was clockwise 
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at all three current meters during four months and counterclockwise during four other 
months. The remaining nine months contained ellipses with opposing rotation. 

The semidiurnal current ellipses at mooring P3 were notably smaller than those 
at P2. The average energy represented by the ellipses at 100 m was 28.4 cm^ sec"^ with 
66% rotating clockwise. At 350 m the average energy decreased to 22.7 cm^ sec^ with 
66% of the currents rotating counterclockwise. A four month average at 500 m resulted 
in an average semidiurnal energy of 34.5 cm^ sec"^ with all four ellipses rotating clock- 
wise. At 1000 m the average energy was only 8.3 cm^ sec"^ At this level, 66% of the 
ellipses rotated clockwise. Like the diurnal currents at P3, the semidiurnal currents were 
greatest at 500 m. The smallest ellipses and least energy was at the 1000 m level. At 
mooring P3 two of the si.x months analyzed had semidiurnal currents at four depths ro- 
tating in a clockwise direction and there was one month when the rotation was 
counterclockwise at each depth. 

The average orientation of 500 m semidiurnal current ellipses at P2 was 030° T 
and the isobaths were crossed at an angle of 25°. The average orientation of 1000 m 
semidiurnal current ellipses at P3 was 043° T and the semi-major axis crossed the 
isobaths at 23°. 

B. ENERGY PROFILES 

Vertical profiles of energy at both the diurnal and semidiurnal frequencies (Figures 
15 and 16) were constructed for each month using the output from the spectral analysis 
program. These plots contain information similar to that of the tidal current ellipses 
expressed in a different format. At P2 the highest energy levels were at the deepest 
current meter for both frequencies. The energy levels at P3 show low energy levels at 
the deeper current meters. 
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!. Spectral Energy Profiles from May 1989 to November 1989 at Mooring 
P2: All plots are scaled identically. 
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Figure 16. Spectral Energy Profiles from December 1989 to September 1990 at 
Moorings P2 and P3: All plots are scaled identically. 
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1. Diurnal Energy Profiles 

At mooring P2 the average monthly diurnal energy levels were 13.9 cm- sec-^ at 
100 m, 13.4 cm^ sec"^ at 350 m and 20.2 cm^ sec"^ at 500 m. The energy at 500 m ex- 
ceeded that at 350 m for each of the seventeen months analyzed except September 1989. 
The energy at 500 m was also higher than that at 100 m except for the months of June 
1989, February 1990 and April 1990. The curves plotted in Figures 15 and 16 exhibit a 
'C' shape with minimum energy levels typically found at a depth of 350 m. 

Diurnal energy at P3 was greatest near the surface. Average monthly values 
were 17.3 cm^ sec^ at 100 m. 7.2 cnF sec^ at 350 m 9.7 cm^ sec"^ at 500 m and 9.8 
cm^ sec“^ at 1000 m. The curves plotted below 100 m are nearly uniform and show vir- 
tually no amplification of energy toward the bottom. 

2. Semidiurnal Energy Profiles 

The semidiurnal current energy levels at P2 were likewise greatest at the deepest 
current meters. The same 'C' shaped curve observed with diurnal currents, with a min- 
imum energy at 350 m, was exhibited by the semidiurnal currents. Average monthly 
values of 63.1, 47.6, and 77.4 cm^ sec^ were found at 100 m, 350 m and 500 m respec- 
tively. The energy at 500 m was greater than that at 350 m for fourteen of the seventeen 
months analyzed. The energy at 500 m exceeded that at 100 m for nine of the seventeen 
months. 

At P3, the average semidiurnal energy levels were fairly consistent until drop- 
ping off sharply at 1000 m. Values of 20.1, 21.3, 30.3 and 7.3 cm^ sec"^ were found at 
100 m, 350 m, 500 m and 1000 m. No bottom amplification of seinidiurnal energy was 
evident at P3. 

C. TEMPORAL VARIABILITY 

The same data used to plot the vertical profiles were used to create time series 
showing the change of energy at each depth for the two moorings. Diurnal bottom 
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trapping was evident at P2, with energy at 500 m greater than the energy at 350 m for 
the entire seventeen month period analyzed (Figure 17). Bottom trapping was always 
present and was most pronounced during July - August 1989, October 1989 - Januart’ 
1990 and May - September 1990 when the energy at 500 m exceeded that of any other 
depth. The magnitude of the diurnal energy at 500 m was statistically greater than 350 
m or 100 m at 95 percent confidence for the period October 1989 through January 1990. 
At mooring P3 no diurnal bottom trapping was indicated by the energy at 350 m, 500 
m or 1000 m (Figure 18). Diurnal energy did not increase toward the bottom and the 
energy near the surface, at 100 m, was always greater than at deeper current meters. 

Bottom trapping of semidiurnal energy occurred at P2 for fourteen of seventeen 
months when the energy at 500 m exceeded the energy at 350 m (Figure 19). In addition, 
the energy at 500 m was greater than 100m for nine of the seventeen months analyzed. 
Bottom trapping at 500 m was statistically significant using a 95 percent confidence in- 
terval for September 1989, December 1989, April - May 1990 and September 1990. At 
P3 no semidiurnal trapping was evident and the energy at 1000 m was less than at any 
other depth for five of the six months analyzed (Figure 20). 

D. APPLICATION OF TRAPPING THEORY 

The size of the current ellipses and energy distribution at each depth suggest that 
an amplification of energy occurs near the bottom of mooring P2. A similar pattern was 
not obvious at P3. To compare observed motions to theoretical trapped waves an ex- 
ponential curve, representing the decay of trapped energy, was compared to the observed 
profile of kinetic energy. The theoretical curves were constructed using equation (12) 
with K given by equation (16). 

1. Trapping Scales and Decay Curves 

At mooring P2 the maximum trapping frequency is exactly equal to the diurnal 
frequency {<o^ = 1 cpd) when, by equation (7), N = 1.930 x 10"'' cps . Using // = 800 m 
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Figure 17. Time Series of Diurnal Energy at Mooring P2: The temporal vari- 
ability of kinetic energy is shown with curves representing the energy 
levels at current meter depths. 
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Figure 18. Time Series of Diurnal Energy at Mooring P3: As indicated in Figure 
17. 
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Figure 19. Time Series of Semidiurnal Energy at Mooring P2: 
Figure 17. 
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Figure 20. Time Series of Semidiurnal Energy at Mooring P3: 
Figure 17. 
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and /= 1.375 X 10"^ cps in equation (16) results in a trapping coefficient of 
K = 3.31 16 X 10"^ nr'. The horizontal length scale of this wave is L = = 26.632 km. 

Similar calculations were conducted for the first five baroclinic modes at both P2 and 
P3 for diurnal and semidiurnal trapping using /V values determined by equation (7). The 
coefficients are tabulated in Tables 3 and 4. 



Table 3. TRAPPING COEFFICIENTS AT MOORING P2 



vj^ = 1 cpd 


iu^ = 2 cpd 


K, = 3.311x10-’ m-‘ 


K, = 1.472x10-’ m-' 


Kj = 6.622x10-’ nr' 


K, = 2.943x10-’ m-' 


Kj = 9.934x10-’ m-' 


Kj = 4.415x10-’ m-' 


1.321x10-’ m-' 


Ki = 5.887x10-’ m*' 


K; = 1.652x10-’ m-' 


K, = 7.358x10-’ m-' 



Table 4. TRAPPING COEFFICIENTS AT MOORING P3 



(i)^ = 1 cpd 


= 2 cpd 


K, = 6.622x10-’ m-' 


K, = 2.949x10-’ m-’ 


K,= 1.324x10-’ m-' 


Kj = 5.887x10-’ m-' 


Kj = 1.986x10-’ m*' 


Kj = 8.830x10-’ m-> 


K 4 = 2.643x10-’ m-' 


K, = 1.177x10-’ m-' 


K< = 3.331x10-’ m-' 


K<= 1.476x10-’ m-' 



For the baroclinic modes, the trapping coelficients in Tables 3 and 4 were used 
to produce a vertical trapping scale. The trapping scale (or e-folding scale) indicates the 
height off the bottom where the energy of a trapped wave decayed by a factor of e~' . 
Most of the trapped energy for each wave was contained in a bottom layer with a 
thickness (determined by the trapping scale) indicated in Table 5. Note that the two 
deepest current meters at P2 (located 450 m and 300 m from the bottom) fell within this 
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layer but the lowest instrument at P3 (located 800 m from the bottom) did not. For 
higher modes trapping is stronger, trapping scales are smaller and energy decays more 
rapidly from the bottom. 



Table 5. TRAPPING SCALES AT MOORINGS P2 AND P3 



Mooring 


= 1 cpd 


o), = 2 cpd 


P2 


302 m 


679 m 


P3 


151 m 


339 m 



Using the appropriate trapping coefficient, exponential decay curves for each 
mooring were constructed which represent the first two baroclinic modes at diurnal and 
semidiurnal frequencies at moorings P2 and P3 (Figures 21 and 22). Similar curves were 
fit to observed spectral energy profiles to determine if the vertical energy distribution off 
Point Sur indicated the presence of bottom trapped waves. The curves w'ere constructed 
using iV values determined from CTD data whenever the diurnal or semidiurnal fre- 
quencies were less than the maximum trapping frequency. 

2. Analysis of Bottom Trapping 

CTD data from each of the hydrographic cruises in Table 2 was analyzed to see 
if bottom trapping was possible at P2 or P3 at diurnal or semidiurnal frequencies using 
the criteria that = N sin a must be equal to or greater than the tidal frequencies. 
Each CTD cast taken at a mooring site was assumed to be representative of the aver- 
aging interval used to compute the spectral estimates (one month). 

When bottom trapping w’as possible, theoretical decay curves were constructed 
using observed values of H, and appropriate trapping coefficients (k). The curves 
were compared to the observed energy levels obtained from spectral analysis by match- 
ing the curves at the level of the deepest current meter (500 m at P2) and comparing the 
theoretical and observed energy at other (100 m and 350 m at P2) levels. When the 
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Figure 21. Theoretical Trapping Profiles at Mooring P2; The exponential decay 
of trapped energy (with arbitrar>’ units) is illustrated for the first and 
second modes at tidal frequencies. 



modeled energy fell within the 95 percent confidence intervals of the observations, the 
change of energy with depth was considered to be consistent with bottom trapped waves. 
The buoyancy frequencies N, maximum trapping frequencies o><, and trapping 

coefficients k, for each month when data were available are tabulated in Table 6. The 

/ 
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Figure 22. Theoretical Trapping Profiles at Mooring P3: The exponential decay 
of trapped energy (with arbitrar>' units) is illustrated for the first and 
second modes at tidal frequencies. 



maximum trapping frequency <u^ at P2 exceeded the diurnal frequency (<u = 1.0 cpd) 
each of the eight months. The maximum trapping frequency exceeded the semidiurnal 
frequency (co = 2.0 cpd) except for the month of January 1990 (co^ = 1.22 cpd). 
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Table 6. OBSERVED BUOYANCY FREQUENCIES, TRAPPING FREQUEN- 
CIES AND MODE 1 TRAPPING COEFFICIENTS AT MOORING P2 



Month 


N (cpd) 


(j), (cpd) 


Kd, ( m-') 


^ StmiJiunial ) 


M.AY 89 


66.2 


3.97 


3.30555 X 10-3 


6.61116 X 10-3 


JUL 89 


43.0 


2.57 


3.30556 X 10-3 


6.61127 X 10-3 


SEP 89 


39.4 


2.36 


3.30557 X 10-3 


6.61130 X 10-3 


NOV 89 


34.7 


2.08 


3.30558 X 10-3 


6.61137 X 10-3 


JAN 90 


20.3 


1.22 


3.30564 X 10-3 


No Trapping 


.MAR 90 


51.8 


3.10 


3.30556 X 10-3 


6.61121 X 10-3 


MAY 90 


38.7 


2.32 


3.30557 X 10-3 


6.61131 X 10-3 


JUN 90 


55.7 


3.34 


3.30556 X 10-3 


6.61120 X 10-3 


AVERAGE 


43,7 


2.62 


3.30557 X 10-3 


6.61126 X 10-3 



The diurnal energy distribution at P2 (Figures 23 and 24) agreed with the mod- 
eled decay for seven of the eight months. During January 1990, the observed energy at 
350 m was less than indicated by theory. The theoretical curves, extrapolated to a depth 
of 800 m, indicated an average kinetic energy level of 40.5 cm^ sec"^ trapped at the sea 
floor. 

The maximum trapping frequency o), at P2 exceeded the semidiurnal frequency 
(o), = 2.0 cpd) during seven of the eight months. The observed energy at 100 m during 
May 1989 was higher than was modeled by matching the theoretical curve to the 500 
m observation. The modeled curve fell within the 95 percent confidence intervals for the 
six remaining months. Extrapolation of the theoretical curve to 800 meters produced 
an average kinetic energy level of 310.0 cm^ sec"^ at the bottom over the time period 
studied. 

The buoyancy frequency calculated over the deepest 10 m of CTD data ranged 
from 20.3 cpd to 66.2 cpd and the magnitude of the maximum trapping frequency, which 
was proportional to N , ranged from 1.22 cpd to 3.97 cpd. The fluctuations of had a 
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much smaller efiect on the trapping coefTicients. The diurnal trapping coefTicient at P2 
varied only 9.0 x 10"® m"‘. The semidiurnal trapping coefTicient had a range of only 
2.1 X 10"'' m"‘. The consistent trapping coefficients produced nearly identical trapping 
scales at P2. The diurnal trapping scale (defined by the e-folding scale) was 302 m and 
the semidiurnal trapping scale was 151 m for each of the eight months analyzed. 

The diurnal energy distribution at P3 did not have a profile indicative of bottom 
trapped waves. Since the maximum trapping scale for a diurnal wave at P3 was ap- 
proximately 151 m and the deepest current meter at P3 was 800 m above the sea floor, 
any energy from a trapped wave would have decayed 99.5 percent before reaching this 
current meter and could not be observed. The maximum trapping frequency at this 
mooring never exceeded (o^= 2.0 cpd and no semidiurnal trapping occurred. 
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line indicating 2.0 cpd in the top row. The middle row shows obser%’ed 
diurnal spectral energy profiles (solid lines) with 95% confidence inter- 
vals (horizontal bars). Theoretical energy profiles (dashed lines), based 
on a decay scale calculated using A' from the bottom of CTD casts (or 
at 800 m for July 1989) , are matched to the observations at 500 m. 
Similar profiles are shown in the third row for semidiurnal energy. 
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plotted as described in Figure 23. No theoretical semidiurnal energy 
profile is shown for Januar>- 1990 when A' sin a near the bottom was 
less than 2.0 cpd. 
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V. DISCUSSION 



a\. bottom trapped waves 

The trapping or amplification of tidal band kinetic energy near the ocean bottom 
was evident from the data obtained at mooring P2. Both the current ellipses and the 
energy distribution displayed a significant increase from a depth of 350 m to 500 m 
during the seventeen months of data analyzed. A similar increase was not evident from 
the data obtained at mooring P3. 

The two methods of analysis used to examine the currents produced results which 
were in general agreement. The energy levels derived through complex demodulation 
and current ellipses were larger in magnitude than those obtained through spectral 
analysis. The total KE represented by the ellipses was 10.7 % greater than the energy 
obtained from spectral analysis. The lower energy spectral analysis values may be due 
to Doppler shifting of small wavelength internal tides to other parts of the spectrum. 
This shifting could spread energy around the tidal band so that energy measured pre- 
cisely at the frequency of a tide, using the high resolution FFT of VCSPC3, improperly 
excluded tidal energy located in sidebands {IVunsh, 1975). The difference between the 
energy from ellipses and spectral analysis was greatest at P3 during December 1989. The 
spectral density plots for the current meter records at 100 m and 350 m show that the 
diurnal peak was spread out over a wider frequency band (perhaps by Doppler shifting) 
and had a smaller amplitude than for other months of the same record. 

Both ellipses and spectral energy profiles indicated more energy was present at 100 
m than at 350 m at both moorings for diurnal and semidiurnal currents. The current 
meters at 100 m were above the main pycnocline and the higher energy levels in the 
surface layer could be the result of wind forcing at or near tidal frequencies. Strong 
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coastal winds with diurnal variations have often been observed along the California 
coast {Beardsley ei al., 1987). The pycnocline could prevent coupling of surface energy 
to deeper layers. Diurnal tidal currents may also be influenced by inertial currents with 
a frequency is 1.1 cpd off Point Sur. The narrow frequency separation between diurnal 
and inenial currents could result in smearing or leakage of energy into the diurnal band. 

The observed distribution of tidal energy had a vertical decay rate indicative of 
bottom trapping at mooring P2. The fit of a theoretical curve to the data fell within the 
95 percent confidence intervals with few exceptions. In January 1990 the observed 
diurnal energy at 350 m was less than predicted by bottom trapped theory using the 
buoyancy frequency calculated from a single CTD cast. If the buoyancy frequency had 
been greater the trapping coefficient, k, would have been larger. A larger k would result 
in a faster decay rate and the theoretical energy at 350 m would have been closer to 
observations. The semidiurnal energy at 100 m during May 1989 was higher than ex- 
pected using the modeled decay curve and these energetic semidiurnal currents may re- 
sult from surface (wind) forcing. Another factor in the vertical semidiurnal energy 
distribution is free internal waves w'hich exist at superinertial frequencies. Free internal 
waves could be generated by tidal forcing across the continental shelf break and could 
then propagate into deeper stratified water {Baines, 1982). Horizontal velocities over the 
shelf have been observed w'here the highest velocities were near the surface and bottom. 
The velocity structure was attributed to the first dynamical mode (Rosenfeld, 1990) and 
the same structure may dominate internal tides over the continental slope. The 
semidiurnal currents at P2, located close to the shelf break, were twice as energetic as the 
semidiurnal currents found further offshore at P3 and did have more energy at the top 
and bottom current meters (100 m and 500 m) than at the middle (350 m). Bottom 
trapping did not occur with any discernable seasonal pattern. The absence of distinct 
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temporal variation agrees with other observations of internal tides along the California 
coast {Noble et ai, 1987). 

The theoretical curves plotted for P2 in Figures 23 and 24 were extrapolated to 800 
m. The average diurnal energy represented by these estimates was 40.5 cm^ sec"^ Using 
the average diurnal current ellipse eccentricity (e = 3.42), diurnal currents at 800 m 
would have had a semi-major axis of 6.1 cm sec~‘ and a semi-minor axis of 
1.8 cm sec"‘. The theoretical energy for semidiurnal currents extrapolated to 800 m was 
310.0 cm^ sec”^ An average semidiurnal eccentricity of c = 3.26 produced an ellipse re- 
presenting bottom currents with semi-major and semi-minor axis lengths of 16.8 
cm sec”' and 5.2 cm sec“‘ respectively. The magnitude of these currents is reasonable, 
however, bottom trapped theory' assumes inviscid boundary' conditions and bottom 
friction may reduce the velocity of predicted bottom trapped currents. 

The average wavelength of bottom trapped diurnal waves, computed using the av- 
erage trapping coefficient in Table 5, was 66.5 km. The average wavelength of bottom 
trapped semidiurnal waves was 33.2 km. This compares favorably to first baroclinic 
mode internal waves of semidiurnal frequency observed during CODE which had wave- 
lengths between 20 and 30 km (Rosenfeld, 1990). 

The rotation of tidal current ellipses is not described by bottom trapped theory'. The 
observed rotation of the current ellipses was variable with 73 percent of the diurnal el- 
lipses rotating counterclockwise and 52 percent of the semidiurnal ellipses rotating 
clockwise. Barotropic Keh in waves, free internal waves and shelf waves all influence the 
rotation of current ellipses at tidal frequencies. The observed direction of rotation is 
determined by the cumulative effect of these unresolved wave regimes. 

Barotropic semidiurnal tides are dominated by the A/j constituent. The horizontal 
velocities generated by Kelvin waves at this frequency are normally aligned with topog- 
raphy {Noble et ai, 1987). The observed tidal current ellipses all had cross- slope com- 
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ponents and vertical velocity structures which may reflect the influence of a baroclinic 
wave structure. At the deepest current meters, the diurnal current ellipses were oriented 
more across the bottom slope {P2 = 44°, Pl> = 54°) than were the semidiurnal ellipses ( 
P2 = 25°, P3 = 23° ). With this geometry' the particle motion of diurnal currents had a 
larger upslope component and more bottom trapping than semidiurnal currents of the 
same magnitude. Even though observed diurnal currents were significantly less energetic 
than semidiurnal currents, diurnal bottom trapping was always evident at mooring P2. 

Both moorings had a similar bottom slope a and the currents analyzed were limited 
to diurnal and semidiurnal frequencies. The only variable in determining the limiting 
trapping frequency ( co, = N sin a ) was the buoyancy frequency. Higher values of N 
resulted in higher trapping frequencies. For each of the eight months when CTD casts 
were made at P2, diurnal trapping was possible as indicated by the calculated values of 
(o^ and semidiurnal trapping was indicated during seven of these months. The current 
ellipses and energy profiles exhibited bottom trapping in each case. No bottom trapping 
was evident at P3 although there were two months with (o^ values high enough to permit 

the trapping of diurnal energy. The absence of any notable bottom trapping at P3 may 

« 

be due to the trapping scale and the depth of the instruments at that mooring. All three 
current meters at P3 were above the e-folding depths of bottom trapped diurnal and 
semidiurnal energy. If any trapped energy was present, the magnitude of the signal was 
indiscernible at the P3 current meters. 

The difference in maximum trapping frequencies at each mooring is caused by dif- 
ferent hydrographic conditions at P2 and P3. The cruise of November 1989 can be used 
to illustrate typical conditions. Density gradients are greater near the ocean bottom at 
P2 than at P3 (Figure 25) and the gradients produce higher N values at P2 than those 
calculated at P3. The stronger gradients at P2 permit trapping ofw’aves with frequencies 
higher than the maximum frequency of trapped waves at mooring P3. 
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Figure 25. Vertical Density Section Off Point Sur, California During November, 
1989: The positions of moorings P2 and P3 are indicated. 

B. LIMITATIONS 

Several assumptions were made in comparing the observations at P2 with theoretical 
predictions for bottom trapped waves. The first was that the tidal currents were always 
oriented across the sea lloor slope so that the angle 4> = 90'’. This angle permits maxi- 
mum trapping and simplifies the equation cr = /V sin a sin 0 to a>, = Vsina. At angles 
less than 90°, the limiting trapping frequency o)^ will be lower and the bottom trapping 
of tidal energy may not occur. If </> was assumed to be the angle at which the semi-major 
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axis of the tidal current ellipses crossed the bottom contours, then the maximum trap- 
ping frequency a at P2 (0 = 44° at 500 m) would still exceed 1.0 cpd for seven of the 
eight months when N values were available. For the P2 semidiurnal ellipses ( 
4> = 25° at 500 m), o would not have exceeded 2.0 cpd for any month. If no semidiurnal 
trapping was possible, then the vertical distribution of energy at semidiurnal frequencies 
must be attributed to other factors including free internal waves. 

A second assumption was that a single CTD cast provided information which could 
be used to interpret monthly averaged current meter data. The CTD casts typically 
ended 50 m above the sea floor and salinity and temperature values over the deepest 10 
m interval were used to calculate N for each month. These N values were assumed to 
be constant from the bottom of the CTD cast to the sea floor. Trapping frequencies 
calculated with N were similarly extrapolated. The use of extrapolated values was sup- 
ported by the CTD cast of July 1989 which was slightly offshore of mooring P2 and ex- 
tended to 903 m (the 10 m interval centered at 800 m was used to calculate N at P2). 
The values of N were fairly constant during July 1989 (42.2 cpd - 40.9cpd) at the depths 
where extrapolation was used (650 m to 800 m) for the seven other CTD casts. As dis- 
cussed in the previous section, variations in iV could affect the fit of the theoretical decay 
curves to observations. The a>^ and k values using a constant N were compared to 
co^ and k computed using N obtained by extending the slope of the curves to 800 m. 
The differences created slight variations of which did not conflict with any of the re- 
sults presented in Chapter IV. The differences in the k values were very small and con- 
sidered insignificant. 

The instrumentation at each of the moorings was not ideally positioned for the study 
of bottom trapped waves. The current meters were too sparse to provide a good profile 
of the currents throughout the water column. The current meters at the 100 m level were 
above the pycnocline and may have contained energy from surface forcing which was 
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not considered in constructing the theoretical decay curves. The lowest current meters 
at P2 and P3 were not deep enough to detect the strongest currents expected near the 
bottom. At P3 the small trapping scales made it unlikely that any signal from a bottom 
trapped wave could be detected. The current meter at 1000 m was well above the e- 
folding depth of diurnal and semidiurnal waves. Additionally, the location of P3 in an 
area with relatively weak density gradients restricted most bottom trapping to subtidal 
frequencies. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The analysis of current meter records from moorings located off Point Sur, 
California showed an increase of energy at tidal frequencies with depth. This can po- 
tentially be explained using the theor>’ of Rhines (1970) for bottom trapped waves over 

> 

the continental slope. The increased energy was observed whenever bottom slope, water 
stratification and frequency met prescribed conditions. 

Seventeen months of current meter data was examined from mooring P2. The data 
were examined using complex demodulation and spectral analysis. At mooring P2, 
bottom trapped energy at the diurnal frequency was observed for each month analyzed. 
Bottom trapping occurred at semidiurnal frequencies for fourteen of the seventeen 
months but the energy contribution due to free internal waves, which influenced the 
semidiurnal results, was not determined. Bottom trapping at mooring P2 occurred dur- 
ing each month of the year and no seasonal patterns could be distinguished. The maxi- 
mum trapping frequency at P2, computed using N sin a, was an accurate indicator of 
bottom trapping at tidal frequencies. 

No bottom trapping was observed at mooring P3 using data collected over six 
months. Weak density gradients near the sea floor often kept the maximum trapped 
frequency below the tidal frequencies. The current meters at this mooring were posi- 
tioned too far above the sea floor for detection of bottom trapped waves at diurnal or 
semidiurnal frequencies. 

The bottom trapped wave theory' described by Rhines (1970) was also used to cal- 
culate vertical decay scales and energy distribution profiles. The theoretical decay was 
modeled using hydrographic data collected near the mooring locations. In most cases 
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the modeled decay agreed with the observed decay of tidal energy away from the sea 
floor within the error of the calculation. 

The tidal currents off Point Sur are the most energetic currents found over the con- 
tinental slope. The semidiurnal currents dominate but diurnal currents can be signif- 
icant, particularly when amplification due to bottom trapping occurs. Further study of 
bottom trapping should be conducted using the continuing current meter arrays at or 
near the location of P2. Additional deep instrumentation, particularly just above the sea 
floor, is needed to adequately study this phenomena. Because trapping was observed 
throughout the year, a short deployment of a few months length may collect sufficient 
data to verify trapping at tidal frequencies. Tidal currents are at the upper limit of the 
range of trapped waves and longer deployments could be used to examine coastally 
trapped waves occurring at much lower frequencies under the same conditions. 
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